Angular Broadening of Intraday Variable AGN. II. Interstellar 

and Intergalactic Scattering 



T. Joseph W. Lazio 

Naval Research Laboratory, 4555 Overlook Ave. SW, Washington, DC 20375-5351; USA 

Joseph. LazioOnrl .navy .mil 

Roopesh Ojha 

NVI/United States Naval Observatory, 3450 Massachusetts Ave. NW, Washington DC 20392-5420; USA 

ro j haSusno . navy . mil 

Alan L. Fey 

United States Naval Observatory, 3450 Massachusetts Ave. NW, Washington DC 20392-5420; USA 

af eySusno . navy .mil 

Lucyna Kedziora-Chudczer 
Institute of Astronomy, School of Physics A28, The University of Sydney, NSW 2006; Australia 

James M. Cordes 

Department of Astronomy, Cornell University and NAIC, Ithaca, NY 14853; USA 

cordesOastro . Cornell . edu 
and 

David L. Jauncey and James E. J. Lovell 

Australia Telescope National Facility CSIRO, PO Box 76, Epping, NSW 1710; Australia 



ABSTRACT 

We analyze a sample of 58 multi-wavelength, Very Long Baseline Array observations of active 
galactic nuclei (AGN) to determine their scattering properties. Approximately 75% of the sample 
consists of AGN that exhibit centimeter- wavelength intraday variability (interstellar scintillation) 
while the other 25% do not show intraday variability. We find that interstellar scattering is 
measurable for most of these AGN, and the typical broadening diameter is 2 mas at 1 GHz. We 
find that the scintillating AGN are typically at lower Galactic latitudes than the non-scintillating 
AGN, consistent with the scenario that intraday variability is a propagation effect from the 
Galactic interstellar medium. The magnitude of the inferred interstellar broadening measured 
toward the scintillating AGN, when scaled to higher frequencies, is comparable to the diameters 
inferred from analyses of the light curves for the more well-known intraday variable sources. 
However, we find no difference in the amount of scattering measured toward the scintillating 
versus non-scintillating AGN. A consistent picture is one in which the scintillation results from 
localized regions ("clumps") distributed throughout the Galactic disk, but which individually 
make little contribution to the angular broadening. Of the 58 AGN observed, 37 (64%) have 
measured redshifts. At best, a marginal trend is found for scintillating (non-scintillating) AGN 
to have smaller (larger) angular diameters at higker redshifts. We also use our observations to try 
to constrain the possibility of intergalactic scattering. While broadly consistent with the scenario 
of a highly turbulent intergalactic medium, our observations do not place significant constraints 
on its properties. 



1. Introduction 

There is now compelling evidence that the in- 
traday variability (IDV) phenomenon — intensity 
variations on hour time scales at centimeter wave- 
lengths in coi npact, flat-spec t rum, extragalactic 



sourc es (e.g., iHeeschenI Il984t iQuirrenbach et al 
19921 ) — is of extrinsic origin. Density fluctua- 
tions in the interstellar medium (ISM) induce 
refractive index fluctuations which, when com- 
bined with the relative motions of the scatter- 
ing medium and the Earth, produce intensity 
variations or scintillations. The evidence for 
this conclusion is both differences in the vari- 
ability pattern arrival times at widely spaced 
radio telescopes and annual cycles in the vari- 



(Jauncev et all 2000 
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20061 Bienall et al.l2006h. 



In order to exhibit interstellar scintillations 
(ISS), a source must contain a sufficiently compact 
component (analogous to "Stars twinkle, planets 
don't" ) such that its angular diameter is compara- 
ble with or smaller than the size of the first Fresnel 
zone of the scattering screen, i.e., of order tens of 
microarcse conds at frequ encies near a few giga- 
hertz (e.g.. lWalke3ll998h . The absence of ISS in 
a source could be either an intrinsic or extrinsic 
effect. Active galactic nuclei (AGN) might have a 
range of intrinsic diameters, in which case only the 
most compact would exhibit ISS. Alternately, in- 
terstellar density fluctuati ons produce a rich range 
of observable phenomena (|Rickettlll990[ ). of which 
scintillations are only one manifestation. Angular 
broadening along the line of sight, due either to 
multiple ionized media or an extended medium, 
could produce apparent diameters of AGN suffi- 
ciently large that the AGN would not display ISS. 

Consistent with this re quirement f o r a co mpact 
(~ 10 fj,as) component, lOiha et al.l (|2004l ) have 
compared AGN that display ISS with those that 
do not and find that the scintillating AGN typ- 
ically are more core-dominated on a milliarcsec- 
ond scale than the non-scintillating AGN. This 
result is striking given that their observations com- 
pared the source structure on milliarcsecond, not 
microarcsecond, scales. However, the observations 



of 8.4 GHz. This frequency is sufficiently high that 
interstellar scattering effects on most lines of sight 
through the ISM would not be detectable on ter- 
restrial baselines nor could they have exploited the 
wavelength dependence for scattering to separate 
the its effects from the intrinsic diameters of the 
sources. 

The AGN observed by lOjha et all ([2004) were 
drawn from the Micro- Arcsecond Scintillation- 



Induced Variability (MASIV) survey (jLovell et al 
2OO3I ). The MASIV survey is a large variability 



survey of the northern sky with the primary goal 
being the construction of a large sample of scin- 
tillating AGN. The survey used the Very Large 
Array (VLA) at 5 GHz in a multi-array mode and 
has yielded scintillation information on over 500 
AGN, of which over half have been found to be 



scintiUating (jLovell et al.ll2007n . 

Ojha et al.l (|2006h presented multi-frequency 
observations of a subset of MASIV sources. Their 
observations were designed to create a sample of 
sufficient size to compare and contrast the scatter- 
ing behavior of scintillating and non-scintillating 
AGN. This paper reports the analysis of those 
observations. 

2. Source Sample 

Our sample consists of 58 AGN, observed in 
two different programs. The first subset con- 
sists of 49 AGN from the MASIV survey, ob- 
served over 3 days with the Very Long Baseline 
Array (VLBA) in 2003 February At the time 
of the observations, approximately half of these 
49 AGN were classified as highly variable MASIV 
sources, with scintillation indices larger than 2%; 
the other half were classified as non-scintillators, 
with no scintillation index larger than 0.2%. These 
sources were chosen without regard to their Galac- 
tic latitude or longitude. Since our observations, 
however, further MASIV observations and analy- 
sis shows that many of the AGN identified origi- 
nally as non-scinti llating are in fact scintillating 



of lOiha et al.l (|2004l ) were at the single frequency 



([Lovell et al.l 120071 ). The number of recognized 
scintillating AGN is now 46, and the number of 
non-scintillating sources is now 13 (ratio of 3:1). 
The subset of AGN from the MASIV survey were 
observed at 0.33, 0.61, and 1.6 GHz, with addi- 
tional observations at 2.3 and 8.4 GHz. In total, 
194 VLBA images of the 49 MASIV extragalac- 
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tic radio sourc es at up to 7 obs erving: frequencies 
were obtained ( Ojha et al. 20061 ). Additional data 
were obtained from the United States Naval Ob- 
servatory (USNO) Radio Reference Frame Image 
Databas^ll (RRFID) and the literature. 

The observations of these 49 AGN were ac- 
quired by cycling through the sources so as to in- 
crease the u-v plane coverage. Typical times on 
source range from 10 min. at the higher frequen- 
cies to 25 min. at the lower frequencies. Typical 
noise levels were within a factor of 2-3 of the ex- 
pected thermal noise limits. 

Ojha et al. ( 20061 ) fit gaussian component mod- 



els to the visibility data of the sources, using the 
images as guides. If more than one component was 
required to model a source at a given frequency, 
the most compact component was identified as the 
"core," as the most compact component will be the 
one for which scattering effects will be most appar- 
ent. The most compact component is frequently 
the brightest one. For the few sources where this 
does not hold strictly at all frequencies, the com- 
pact and bright component that could consistently 
be identfied as the same at all frequencies was 
identified as the "core," e.g., J0713-h4349 where 
the northernmost component is identified as the 
core even thoug h it is not the br ightest compo- 
nent at 8.4 GHz (lOjha et al.ll2006l ). 

The second subset consists of 9 AGN cho- 
sen from those used to de fine the International 
Celestial Reference Frame ( Johnston et al. 19951 : 



Ma et al.lll998f ). The initial motivation was to use 



these AGN to search for scattering resulting from 
the interstellar media of galaxies along their lines 
of sight. The sources were chosen to be (1) At 
Galactic latitudes |6| > 45°; (2) Strong, with 
Seem ~ 1 Jy; and (3) Compact and domi nated 
by a single component (jFev fc Charlotlll997l l. Ob- 
servations of these 9 sources were conducted with 
the VLBA on 2001 February 17 and 18 at 0.33 
and 0.61 GHz. Calibration, imaging, and model 
extraction w as performed in a manner similar to 
that used bv lOiha et all (|200d ). 

Because both subsets involved observations 
at 0.33 and 0.61 GHz, the observations were gen- 
erally carried out at night so that the sources 
had large solar elongations. Indeed, observing at 
large solar elongation was an explicit criterion in 



scheduling the observations for the second subset. 
As a result, the smallest elongation for any source 
is 75°, and the typical elongation is approximately 
130°. 

For the present analysis, we have used the core 
components models from these two observing pro- 
grams, augmented by measurements from the lit- 
erature. All sources have angular diameters mea- 
sured at at least 3 frequencies, and some sources 
have measured angular diame ters at as many a s 
7 frequencies. See Table 2 of lOiha et all (|2006l ). 
While the selection criteria for the two subsets dif- 
fered, the sources were treated identically in the 
following analysis. 

3. Analysis 

Angular broadening is manifested as an ob- 
served angular diameter that scales approximately 
as A^. We fit the measured angular diameters to 
the functional form 



(1) 



^ http : //www . usno . navy . mil/RRFID/ 



where 9s and 9i are the scattering and intrinsic 
(FWHM) diameters of the AGN, respectively, at 
the fiducial frequency of 1 GHz. We found the 
best-fitting values for Og and 9i in a minimum 
sense. We considered both a; = (i.e., frequency- 
independent intrinsic diameter, for a flat spectrum 
source) and x = —1 (i.e., frequency scaling for 
a single incoherent synchrotron component) and 
selected the value of x that produced the lower 

As a motivation for the use of equation ([T]), 
as well as anticipating later discussion, we be- 
gin by considering a crude approximation to equa- 
tion ID) . For the sources having measured angular 
diameters at both 0.33 and 1.6 GHz, we have as- 
sumed a simple power-law scaling for the angular 
diameter, 6 (x i/^, and solved for f3. We chose 
0.33 GHz because the frequency dependence of 
scattering means that it will be the most sensi- 
tive to scattering; we chose 1.6 GHz as the second 
frequency as an attempt to balance between hav- 
ing a sufficiently large frequency dynamic range so 
as to obtain a robust estimate of f3 but not having 
such a large frequency range that intrinsic struc- 
ture might dominate. If scattering is important, 
we expect /3 ~ — 2. We find an average value of 
(3 — —1.9 ± 0.1 for the entire sample. Clearly, we 
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anticipate that intrinsic structure may be impor- 
tant for some sources, but that the P is close to 
the expected value from scattering indicates that 
scattering is important for the sample of sources 
as a whole. 

Table [1] summarizes the inferred scattering and 
intrinsic diameters from fitting the data for each 
source to equation |T]). For comparison, we also 
list the predicte d diameter at 1 GHz from the 
NE2001 model (|Cordes fc Laziol |2006|) for inter- 
stellar scattering, 6'ne20oi- Figure [T] illustrates ex- 
amples of the measured angular diameters, show- 
ing the results for an AGN for which a relatively 
large amount of scattering is inferred and one for 
which a relatively small amount of scattering is 
inferred. 
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Table 1 

Fitted Scattering and Intrinsic Diameters 



Name 


i 

n 


b 

n 


z 


Scintillate? 


X 


(mas) 


(mas) 


N 




^NB2001 

(mas) 


J0102+5824 


124.419 


-4.435 


0.644 


Y 


1 


3.6 


1.2 


6 


6 


4.9 


J0217+7349 


128.927 


11.964 


2.367 


N 


1 


2.8 


0.5 


9 


11 


2.0 


J0343+3622 


157.530 


-14.691 


1.484 


Y 


1 


3.2 


3.7 


5 


39 


1.9 


J0349+4609 


152.152 


-6.369 




N 


1 


2.1 


1.2 


4 


5 


4.1 


J0403+2600 


168.025 


-19.648 


2.109 


Yi 





2.4 


0.3 


8 


17 


1.5 


J0419+3955 


160.461 


-7.336 




Y 





2.6 


0.2 


3 


2 


3.8 


J0423+4150 


159.705 


-5.381 


2.277 


Y 





4.1 


0.5 


3 


0.3 


4.6 


J0451+5935 


149.323 


9.660 




Y 





4.9 


0.2 


3 


8 


2.9 


J0502+1338 


187.414 


-16.745 




Y 





3.7 


0.1 


3 


1 


1.6 


J0503+0203 


197.911 


-22.815 




N 





4.6 


0.5 


8 


10 


1.2 


J0507+4645 


161.025 


3.716 




Yi 





1.6 


1.0 


3 


9 


5.2 


J0509+0541 


195.405 


-19.635 




Y 





3.0 


0.2 


5 


24 


1.3 


J0539+1433'' 


191.597 


-8.660 


2.69 


Y 


1 


3.4 


2.4 


3 


3 


3.1 


J0539+1433'' 


191.597 


-8.660 


2.69 


Y 





6.1 


0.4 


3 


3 


3.1 


J0607+6720 


146.804 


20.858 


1.97 


Yi 





1.4 


0.4 


11 


20 


1.3 


J0650+6001 


155.842 


23.155 


0.455 


Y 





7.8 


0.3 


6 


36 


1.3 


J0654+5042 


165.680 


21.106 




Yi 





1.8 


0.2 


3 


45 


1.4 


J0713+4349 


173.792 


22.199 


0.518 


N 


1 


1.3 


4.7 


11 


18 


1.4 


J0721+7120 


143.981 


28.017 


2.06 


Yi 


1 


0.6 


0.9 


8 


39 


1.1 


J0725+1425 


203.643 


13.908 




Y 





1.5 


0.2 


6 


36 


1.7 


J0738+1742 


201.846 


18.070 


0.424 


Yi 





1.4 


0.2 


9 


38 


1.5 


J0745+1011 


209.796 


16.592 


2.624 


Yi 


1 


0.6 


2.2 


10 


36 


1.4 


J0757+0956 


211.311 


19.057 


0.266 


Y 





1.1 


0.2 


9 


24 


1.3 


J0808+4950 


169.163 


32.564 


1.418 


Yi 


1 


1.1 


1.1 


12 


35 


1.1 


J0830+2410 


200.021 


31.876 


0.939 


Y 


1 


0.8 


1.9 


6 


9 


1.2 


J0831+0429 


220.693 


24.331 


0.180 


Yi 


1 


4.8 


1.0 


9 


18 


1.1 


J0842+1835 


207.275 


32.480 


1.270 


Yi 





4.1 


0.5 


8 


23 


1.1 


J0914+0245 


228.352 


32.819 


0.427 


Y 





2.8 


0.1 


4 


16 


1.0 


J0920+4441 


175.700 


44.815 


2.190 


N 


1 


1.7 


0.7 


10 


36 


1.0 


J0956+2515 


205.511 


50.981 


0.712 


Y 


1 


0.0 


2.3 


11 


31 


0.9 


J0958+4725 


170.055 


50.730 


1.882 


Y 





2.3 


0.2 


7 


7 


0.9 


J1008+0621 


233.521 


46.012 




Y 





5.0 


0.0 


3 


8 


0.9 


J1014+2301 


210.699 


54.431 


0.565 


Y 





3.5 


0.3 


7 


37 


0.9 


J1041+5233 


157.521 


54.965 


0.677 


Y 





1.4 


0.8 


4 


1 


0.8 


J1125+2610 


210.920 


70.885 


2.341 


Yi 


1 


0.6 


1.1 


5 


18 


0.6 


J1153+8058 


125.719 


35.836 


1.250 


Yi 





1.9 


0.2 


6 


17 


0.9 


J1159+2914 


199.413 


78.374 


0.729 


Y 


1 


2.7 


1.5 


8 


18 


0.5 



5 




Fig. 1. — Two examples of fitting results. For 
both plots, solid circles show measured diameters, 
arrows indicate upper limits, the solid line indi- 
cates the fit of equation ([1]) to the observations, 
the dashed line indicates the inferred scattering 
diameter, and dotted line indicates the inferred 
intrinsic diameter. Uncertainties on the angular 
diameters are plotted, but in many cases are com- 
parable to the size of the symbol. Top The source 
J2022+6136 for which a relatively large scatter- 
ing diameter, 3.1 mas at 1 GHz, is inferred. Bot- 
tom The source J0745+1011 for which a relatively 
small scattering diameter, 0.6 mas at 1 GHz, is 
inferred. 

There are values of the reduced in Table [1] 
that are larger than unity, at times by a signifi- 
cant factor. These result from a combination of 
two factors. First, the uncertainties for some data 
are likely to be underestimated. For the observed 
diameters obtained from our observations, we es- 
timated their uncertainties using a bootstrap-like 
procedure in which visibility data associated with 
different antennas were removed before perform- 
ing the fit. The range of fitted diameters sug- 
gests a 10% uncertainty. For angular diameters 
obtained from the literature, we have assumed the 
same fractional uncertainty. However, we have 
been able to identify data for which this assumed 



10% is likely to be too small. 

The fitting procedure also yields an estimate 
of the uncertainty in the inferred scattering diam- 
eter (at 1 GHz). The median value of this un- 
certainty is 0.1 mas. We have repeated the fit- 
ting procedure with larger uncertainties adopted 
for the angular diameters obtained from the liter- 
ature, and in some cases even removing apparent 
outliers. The typical change in the inferred diame- 
ter is comparable to the uncertainty in the inferred 
scattering diameter. We have also repeated anal- 
yses described below with larger uncertainties for 
the angular diameters from the literature and find 
no change from the results we present below. 

A second potential cause of large values is 
that we have adopted fixed frequency scaling ex- 
ponents in equation ([T]). We might obtain a bet- 
ter fit by letting a; be a fitted parameter, fitting 
the scattering frequency dependence rather than 
adopting —2.2, or both. To do so would often re- 
quire more data than are available. Consequently, 
while larger than unity, we consider there to be 
plausible explanations for the values in Tabled] 
and shall use the angular diameters resulting from 
our fits. 

divided the sources into 



Oiha et al 



two groups, "scintillators" and "non-scintillators." 
Since the publication of that paper, it has been 
realized that some of the AGN identified as scin- 
tillators displayed variation at only a single epoch 
("once-er") among the MASIV observations, lead- 
ing to the possibility that a non-scintillator would 
be classified mistakenly as a scintillator. Analy- 
sis of the MASIV observations (to be published 
elsewhere) suggests that the light curve from an 
individual epoch can be classified correctly at the 
95% confidence level. From the four epochs of 
observations comprising MASIV, the probability 
of a false identification is only 4.3%, meaning that 
we expect only 1 source in our sample to be clas- 
sified mistakenly. While we identify the single- 
epoch variable sources in Table [H their presence 
should have have no significant effect on our anal- 
ysis, and we treat the single-epoch variable AGN 
as scintillators. Also, analysis has continued on 
the MASIV sources, so there may be occasional 
differences in the classification (scintillating vs. 
non-scintil lating) in our Tab le [T] as compared to 
Table 1 of lOiha et all (|20n()h . 

Figure [2] shows the distribution of the in- 
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Table 1 — Continued 



Name 


e 


b 


z 


Scintillate? 


X 




B^ 


iV 




SnE2001 




n 


n 








(mas) 


(mas) 






(mas) 




O.ooU 


on 


1 /inn 


IN 


1 
i 


l.D 


1 9 


A 


9S 
Zo 


n (\ 

U.D 






i o.ZOi 


U.U / D 


N 
i>i 


1 


n n 

U.U 


o.u 


o 
o 


Q 


n ^ 

U.O 


J1642+6856 


100.705 


36.621 


0.751 


Y 


1 


1.2 


1.7 


13 


46 


0.8 


J IDOD-rDUlZ 


oy .DZ ( 


^7 


n R9^ 

U.DZO 


V 


n 
u 


D.O 


U.O 


Q 
O 


U. 1 


U.O 


J 1 1 4D-i-DZzD 


01 A91 
yl.DZl 


ol.OZU 


o.ooy 


V 

1 1 


A 

U 


9 7 
Z. / 


n 9 
u.z 





1 Q 

ly 


U.O 


J lolZ-t-ODUo 


o4.0o ( 


97 /I7Q 




V 


A 

u 


O.O 


U.O 


Q 
O 


1 9 

iz 


U.O 


J1823+6857 


99.210 


27.669 




Y 


1 


1.8 


1.6 


4 


12 


0.9 


J1927+6117 


92.726 


19.446 




Yi 


1 


1.8 


1.5 


5 


1 


1.0 


J2002+4725 


82.219 


8.793 




Yi 


1 


0.0 


2.7 


3 


4 


1.8 


J2009+7229 


105.355 


20.180 




Y 





1.3 


0.3 


3 


21 


1.1 


J2022+6136 


96.082 


13.775 


0.227 


N 





3.1 


0.4 


15 


30 


1.4 


J2230+6946 


111.248 


10.164 




Y 





2.4 


0.3 


7 


16 


2.0 


J2311+4543 


105.315 


-13.703 


1.447 


Y 





2.4 


0.4 


5 


26 


1.4 


B0955+476 


170.055 


50.730 


1.873 


Y 


1 


2.1 


0.1 


7 


34 


1.0 


B1130+009 


264.364 


57.582 




N 





1.8 


0.9 


4 


51 


0.9 


B1226+373 


147.142 


78.938 


1.515 


N 


1 


2.8 


0.1 


4 


15 


0.5 


B1236+077 


294.112 


70.170 


0.400 


N 


1 


1.4 


0.3 


3 


48 


0.8 


B 1402+044 


343.669 


61.169 


3.211 


N 


1 


2.2 


0.3 


6 


48 


0.6 


B1432+200 


21.387 


65.299 




Yi 




2.3 


0.0 


3 




0.5 


B1459+480 


81.122 


57.419 




Y 





1.6 


1.1 


4 


110 


0.8 


B1502+036 


2.226 


50.254 


0.413 


Yi 




1.0 


2.9 


3 


0.2 


0.7 


B1502+106 


11.381 


54.580 


1.833 


N 


1 


2.4 


0.2 


6 


12 


0.6 



''Both V ^ and v^^ dependences for the intrinsic diameter produced equal ^ . 

Note. — Sources are indicated to be either scintillators (Y), non-scintillators (N), or sources observed to vary 
once and presumed to be scintillators (Yi); x is the spectral index for the frequency dependence of the intrinsic 
diameter, equation (IT|; Qs and Qi are the inferred scattering and intrinsic diameters, respectively; A'^ and are 
the number of data and chi-square in the fit for equation (IT|; and 6nb20oi is the predicted interstellar scattering 
diameter from the NE2001 model. 
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ferred scattering diameters plotted as a function 
of Galactic coordinates. We have considered the 
distribution of the sources on the sky as a function 
of both Galactic latitude b and ecliptic latitude /3 
and the distribution as a function of solar elonga- 
tion. There is no statistically significant correla- 
tion of the inferred scattering diameter with either 
coordinate (strictly, the absolute value of the coor- 
dinate), nor with cos(|5|) or cos(|/3|). (Typical cor- 
relation coefficients are approximately 0.1.) The 
use of the cosine of the angle [cos(|6|) or cos(|/?|)] 
attempts to compensate for the increased amount 
of sky at low latitudes. There is also no corre- 
lation between inferred scattering diameter and 
solar elongation. 

The lack of a correlation of the scattering di- 
ameter with Galactic latitude would appear to be 
at odds with previous determinations that inter- 
stellar scatteringj££reasesj;a£i^^^ latitudes 



(e.g.. iRao fc AnanthakrishnanI 119841 ). However 



while no longitude selection criterion was applied 
in constructing our sample, our sample is nonethe- 
less weighted strongly toward the outer Galaxy. 
We have no sources with longitudes in the range 
—60° < £ < 60°, and only a few at low latitudes in 
the range -120° < £ < 120°. Thus, we attribute 
the apparent lack of a correlation between scatter- 
ing diameter and Galactic latitude as a result of 
having few, essentially no, lines of sight into the 
inner Galaxy. 

While there is no correlation of scattering di- 
ameter with Galactic latitude for the entire set 
of sources, scintillating AGN are consistently at 
lower Galactic latitudes than the non-scintillating 
AGN. Table[2]shows that the average absolute val- 
ues for the Galactic latitudes of scintillating and 
non-scintillating sources differ by nearly 20°. No 
such difference is found for the average (absolute) 
ecliptic latitude. 

Figure [3] shows a histogram of the inferred scat- 
tering diameters. We have used a Kolmogorov- 
Smirnov test to assess whether scintillating AGN 
have a different distribution of inferred scatter- 
ing diameters as compared to the non-scintillating 
AGN. We find no statistical difference: the scat- 
tering diameters of scintillating AGN do not dif- 
fer appreciably from those of the non-scintillating 
AGN. 

Examination of the scattering diameters in- 
ferred for individual objects indicates that some of 




Fig. 2. — The distribution of sources observed as 
a function of Galactic coordinates. The Galactic 
anticenter is at the center of the plot, and longi- 
tude increases to the left. Stars indicate sources 
that scintillate; squares indicate non-scintillating 
sources. The size of the symbol is qualitatively 
proportional to the inferred scattering diameter. 




2 4 6 

Scattering Diameter (mas) 



Fig. 3. — The histograms of scattering diameters. 
The open histogram shows the distribution for the 
scintillating sources; the hatched histogram shows 
the distribution for the non-scintillating sources. 
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Table 2 
Statistical Measures 



Scintillators Non-scintillators 



Galactic latitude (absolute value) 
mean 29?5 ± 3?0 46?3 ± 7?0 

ecliptic latitude (absolute value) 
mean 34?9 ± 4?3 28?8 ± 4?9 

Scattering Diameters (at 1 GHz) 
mean (mas) 2.1 ± 0.2 2.1 ± 0.2 

median (mas) 1.9 2.0 



Redshifts 

mean (mas) 1.38 ± 0.18 1.29 ± 0.31 

median (mas) 1.27 1.40 



the largest scattering diameters result from AGN 
for which no measurements exist below 1 GHz. In 
order that these not bias our result, we removed 
these and repeated the K-S test analysis. There 
is no change in the result, that the scattering di- 
ameters for the scintillating and non-scintillating 
AGN are consistent with having been drawn from 
the same distribution. Both the mean and the 
median scattering diameter for scintillating and 
non-scintillating sources is approximately 2 mas 
(Table El). 

From the entire sample, 37 AGN (64%) have 
measured redshifts. There appears to be little dif- 
ference in the redshift distribution of the scintillat- 
ing and non-scintillating AGN, with the two pop- 
ulations having similar means and medians (Ta- 
ble H]). Figure S] shows the distribution of the 
scattering diameters as a function of redshift. We 
have determined the correlation between the scat- 
tering diameters and redshifts for the entire sam- 
ple, as well as splitting it into the two popula- 
tions, scintillating and non-scintillating. There 
is no correlation of the scattering diameter with 
redshift for the entire sample. There may be a 
marginal correlation, at the 5% confidence level, 
between the scattering diameters and redshift, in 
the opposite sense for the scintillating and non- 
scintillating sources. That is, the scattering diam- 
eters of scintillating (non-scintillating) AGN may 



12 3 4 

Redshift 

Fig. 4. — The distribution of scattering diameters 
as a function of redshift. Stars indicate sources 
that scintillate; squares indicate non-scintillating 
sources. 

become smaller (larger) at higher redshifts. 

4. Discussion and Conclusions 

In one sense, our results are broadly consistent 
with what is known about intraday variability and 
interstellar scattering. In our sample, scintillating 
AGN lie typically at lower Galactic latitudes, con- 
sistent with the notion that the scintillation re- 
sponsible for IDV results from the Galactic ISM 

(CD- 
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In addition to the ISM, other ionized media 
along the hnes of sight to these sources are the in- 
terplanetary medium (IPM) and the intergalactic 
medium (IGM). It is possible that the IPM could 
contribute angular broadening at a level sufficient 
to be detectable in this analysis, particularly given 
the use of 0.33 and 0.61 GHz observations. The 
lack of a correlation with either the ecliptic lati- 
tude (Table [2|) or solar elongation, however, im- 
plies that the IPM makes no detectable contribu- 
tion in these observations. The lack of a strong 
correlation with redshift and the systematically 
lower Galactic latitude of the broadened sources 
suggests that the IGM makes no detectable con- 
tribution to our observations, either. 

In a further effort to differentiate the Galactic 
contribution of scattering from any possible inter- 
galactic contribution, we have searched for pulsars 
within 1° of the AGN in our sample. We find no 
pulsars this close to any of our sources. Given the 
relatively low density of pulsars on the sky, a sig- 
nificantly larger sample of AGN would be required 
in order to make such a comparison. 

On the face of it, the result that scintillating 
AGN are broadened at levels comparable to those 
of the non-scintillating AGN appears to contradict 
the requirement (SJT]) that in order to display scin- 
tillation, a source must be sufficiently compact. 
One possibility is that scattering is, in fact, not 
important at all and that we have mistakenly at- 
tributed the effects of intrinsic structure to angular 
broadening. We regard this as unlikely, given that 
the angular diameters of the sample of sources, 
as a whole, scale approximately as expected from 
interstellar scattering. 

In fact, the magnitude of the estimated broad- 
ening does not appear to be problematic from the 
standpoint of quenching the scintillations. Tak- 
ing 2 mas at 1 GHz as a characteristic scatter- 
ing diameter (Table [2]) , the implied scattering di- 
ameter at 5 GHz is 80 /las, comparable to the 
value derived for a numbe r of the well-known 
scintillating sources (e.g., Rickett et al. 19951: 



the scintillation of all scintillating AGN, we would 
not expect a difference between the average Galac- 
tic latitude of scintillating and non-scintillating 
sources (Table [5]) . That such a difference exists 
suggests that the scintillation for most scintillat- 
ing AGN results from scattering material associ- 
ated with the Galactic ISM. 

A consistent explanation for these results is 
obtained if the scintillation is produced from 
small "clumps" of scattering material, distributed 
thro ughout the Galactic di sk . F or instance, 
IPe^nett-Thorpe fc de Bruvnl (|2003l ) determine 
that a "thin screen" is responsible for the ex- 
treme IDV of J1819-I-3845, with the screen being 
about 10 pc distant and having an internal level of 
scattering measured by = 0.5 m^^"/'^. They do 
not provide a quantitative estimate of the thick- 
ness of this screen, but, in order that the screen be 
"thin," it must surely be the case that its thickness 
is AL < 1 pc. 

Suppose we require that the scattering con- 
tributed by such a clump not make a significant 
contribution to the angular broadening. For illus- 
tration purposes, we adopt the amount of broaden- 
ing contributed by the clump to be 6s.c\ ^ 0.2 mas, 
which would be only a 10% contribution to the 
typical broadening diameter that we measure. The 
resul ting scattering measure SM ( Cordes &: Laziol 
20061 ) is then SM^ ~ lO'^-^ kpc m-20/3. If the 
clump has = 0.5 m~^°/'^, the implied thickness 
is 0.05 pc (^ 10^ AU), which would c ertainly qual 
i fy as "thi n." Mo r eover , following iRickett et al 



(|l995l ) and lRickettI (|2002h . it can be shown that a 
more distant scattering clump tends to produce a 
lower scintillation index. 

Thus, our picture is one in which the Galactic 
disk contains (small) "clumps" of scattering ma- 
terial. Lines of sight through the disk are scat- 
tered by the broadly distributed ionized interstel- 
lar medium, so that AGN over the range of lat- 
itudes that we observe have measurable scatter- 
ing diameters. Some (many?) lines of sight pass 
through one of t hese clumps, and AGN having 



Rickett. Kedziora-Chudczer. fc Jauncevll2002 ; Dennett-ThoppMplcctleefimu^ components are then observed 

to scintillate. However, the clumps are small 



20031 : iBignall et al.ll2003h . 

A notable feature of the most extreme IDV 
sources is that detailed analyses of their light 
curves suggest that the scattering medium respon- 
sible for the scintillation lies quite close to the Sun 
(w 25 pc). If local material was responsible for 



enough that they produce effectively no additional 
broadening. This scenario is also broadly consis- 
tent with the notion of "clumps" of material pro- 



ducin g extreme scattering events (jFiedler et al 



19871 ) and parabolic arcs in pulsar dynamic spec- 
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tra (jHill et all 120051) . 

We can also use the difference between the scin- 
tiUating and non-scintillating sources to set quan- 
titative limits on the amount of radio-wave scat- 
tering contributed by the IGM. We adopt 0.5 mas 
at 1 GHz (sa 3a from Table[2]) as the upper limit on 
the difference in the amount of scattering between 
the two populations. The implied scatt ering mea- 
sure is SM < 10"'' kpc m^^o/s (C ordes fc Laziol 



2006f ). In turn, the scattering measure is given by 



SM = CsmFuID, 



(2) 



where D is the distance, is a fluctuation param- 
eter encapsulating aspects of the microphysics of 
the plasma, Ue is the electron density, and Csm = 
1.8 m"^"/^ cm^ is a constant. For a characteris- 
tic redshift of approximately unity (Figure |4]), the 
equivalent (angular-size) distance is I? « 1.5 Gpc, 
implying Fnl < 10~^°'^ cm~^. 

For a baryonic m atter density Q.f,h^ — 0.127 



( Spergel et al. 20061 ). the intergalactic electron 
density can be no larger than < 2.2 x 



10~^ cm~'^, assuming that helium is f ully ion- 
ized ( Sokasian. Abel, fc Hernquist 20021 ). Thus, 
we require F < 10'^, so as not to violate the in- 
ferred limits on scattering. For reference, in the 
diffuse Galactic ISM, F « 0.2, and in the Galactic 
spiral arms, F 10. In turn, the F parameter is 



(3) 



where C is the normalized second moment of the 
fluctuations, e is the fractional variance in rig 
within the plasma, rj is the filling factor, and £q is 
the largest scale on which the density fluctuations 
occur (or outer scale, if the plasma is turbulent), 
in parsec units. Assuming that C ~ e ~ 1, we 
conclude that Ty^g^^ > 10"^. 

The IGM is tho ught to be permeated by shocks 
(jPave et al.l 120011 ) , which might be expected to 
drive 77 — > 1. Given the larger scales available in 
the IGM, £0^1 Mpc would not be unreasonable. 
We are forced to conclude that the current limits 
on intergalactic scattering, while broadly consis- 
tent with the current picture of the IGM, do not 
yet place significant constraints on its properties. 

While we find no indications of intergalactic 
scattering, future observations are warranted. In 
particular, if a scintillating AGN is found close 



to the line of sight to a pulsar, a comparison be- 
tween the two lines of sight would provide strong 
constraints on the amount of Galactic vs. inter- 
galactic scattering. Also, higher-sensitivity obser- 
vations (e.g., with the very long baseline High Sen- 
sitivity Array or HSA) targeting scintillating AGN 
with larger diameters may provide additional con- 
straints. Many of the AGN with the largest di- 
ameters are not detected at the lower frequencies, 
frequencies at which the VLBA alone has a rel- 
atively low sensitivity. The HSA could be used 
to verify whether these AGN do indeed have such 
large scattering diameters or assess to what extent 
intrinsic structure contaminates the scattering di- 
ameter estimates. 

We summarize our findings as follows. In our 
sample of 58 AGN, approximately 75% of the sam- 
ple exhibit intraday variability (interstellar scintil- 
lation) with the other 25% not showing intraday 
variability. Interstellar scattering is measurable 
for most of these AGN, and the typical broad- 
ening diameter is 2 mas. Scintillating AGN are 
typically at lower Galactic latitudes than the non- 
scintillating AGN, consistent with the scenario 
that intraday variability is a propagation effect 
from the Galactic interstellar medium. The mag- 
nitude of the inferred interstellar broadening mea- 
sured toward the scintillating AGN, when scaled 
to higher frequencies, is comparable to that de- 
termined from analyses of the light curves for 
the more well-known intraday variable sources. 
However, we find no difference in the amount of 
scattering measured toward the scintillating ver- 
sus non-scintillating AGN. A consistent picture 
is one in which the scintillation results from lo- 
calized regions ( "clumps" ) distributed throughout 
the Galactic disk, but which individually make 
little contribution to the angular broadening. In 
our sample, 63% of the AGN have measured red- 
shifts. At best, a marginal trend is found for scin- 
tillating (non-scintillating) AGN to have smaller 
(larger) angular diameters at higher redshifts. Fi- 
nally, while broadly consistent with the scenario of 
a highly turbulent intergalactic medium, our ob- 
servations do not place significant constraints on 
its properties. 
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